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Introduction 


The  earth's  gravitational  potential  can  be  expressed  in  a  spherical  harmonic 
series  where  the  potential  coefficients  are  where  the  overbar  indicates 

full  normalisation  (Kaula,  1966,  Heiskancn  and  Moritz,  1967),  The  potential 
degree  variances  arc  computed  from: 

I 

^  (1) 

■  ::'0 

Estimates  of  the  variation  of  were  obtained  from  gravimetry  by  Kaula  (1966) 
who  suggested  the  following  behavior: 


Another  term  of  interest  is  the  anomaly  degree  varuince: 

I 

(Cf,  +  s/,)  (3) 

■  =' o 

where  y  is  an  average  value  of  gravity  over  the  earth. 

Estimates  of  and  c^  arc  needed  for  several  purposes.  We  can  gain 
information  on  the  density  structure  of  the  earth  by  considering  various  mass 
distributions  that  would  imply  a  potential  coefficient  variation  of  the  form  as 
equation  (2)  (Kaula,  1977).  We  can  also  compute  covariance  functions  of 
quantities  related  to  the  earth's  gravity  field  if  wc  have  models  of  the  anomaly 
degree  variances. 

The  need  for  some  simple  formulas  to  describe  cr^  or  c^  have  led  to  a 
number  of  studies  that  have  fitted  empirically  determined  data  to  various  models. 
Such  procedures  were  followed  by  Kapp  (1973)  who  used  limited  satellite  implied 
information  and  limited  terrestrial  gravity  information.  A  more  refined  model 
was  described  in  Tseheming  and  liapp  (1974)  that  became  widely  used  as  the 
basis  for  covariance  computations  needed  for  the  prediction  of  quantities  dependent 
on  the  gravity  field  using  least  squares  collocation.  Moritz  (1977)  suggested  an 
improved  model  for  the  anomaly  degree  variances  that  was  investigated  by 
Jekeli  (1978)  and  Hein  and  jochemezyk  (1979).  Wagner  and  Colombo  (1978) 
analyzed  altimeter  data  to  imply  some  simple  rules  for  Op  that  depended  on  the 
range  of  I  being  investigated.  Heller  and  Jordan  (1979)  used  the  superposition 
of  white  noise  spherical  shells  at  various  depths  with  various  masses  to  generate 
a  disturbing  potential. 

Recently,  several  questions  have  arisen  that  can  be  answered  with  new  data 
that  is  now  available.  Specifically,  we  are  interested  to  see  if  the  real  world 
potential  coefficients  do  foil  off  as  1/ or  could  it  be  1/f  or  some  other 
power.  Are  the  parameters  of  the  c  model  described  by  Tseheming  and  Rapp 
(1974)  substantially  valid  today?  Is  the  model  suggested  by  Moritz  (1977)  abetter 
fit  to  existing  daUi?  Do  the  suggestions  of  Wagner  and  Colombo  (1978)  and  Heller 
and  Jordan  (1979)  fit  our  new  data? 
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The  data  that  we  will  consider  are  the  potential  coefficients  to  degree  36 
from  the  GEM  lOB  field  (Lerchet.  al. ,  1978),  and  the  potential  coefficients  to 
degree  180  based  on  an  adjusted  set  of  64800  l^x  1°  mean  free  air  anomalies 
described  by  Rapp  (1978). 


Potential  Coefficient  Modelling 

The  first  modelling  will  be  done  with  the  following  equation: 

(4) 

X/ 

so  that  for  Kaula's  rule  A  =  10  and  B  =  2.  Shown  in  Figure  1  are  the  values  of 
a  from  GEM  lOB,  from  Kaula's  rule  of  thumb  and  as  implied  by  the  c^  model 
of  Tscherning  and  Rapp.  These  two  models  have  clearly  more  power  in  them 
than  in  the  actual  gravitational  field. 

We  then  carried  out  least  squares  fits  to  (4)  solving  for  A  and  B  that  best 
fit  the  values  from  GICM  lOB  and  the  180, 180  solution  in  tw'o  different  minimization 
procedures.  The  first  procedure  minimized  the  sum  of  the  squares  of  the  residuals 
while  the  second  procedure  effectively  minimized  the  sum  of  the  squares  of  the 
percentage  residuals.  (A  percent  residual  would  be  100  (adjusted  value  -  observed 
value)/(observed  value).)  This  latter  procedure  is  more  effective  when  fitting  to 
data  having  wide  variations  in  magnitudes.  We  give  in  Table  1  values  of  A  and  B 
for  different  weighting  cases.  We  see  that  the  use  of  weights  as  1/obs^  yields 
poort'r  residual  fits  but  slightly  better  percentage  fits.  The  differences  between 
the  GEM  lOB  and  the  180,180  solutions  do  not  appear  significant.  Plots  of  GEM  lOB 
and  the  two  GEM  lOB  fits  arc  shown  in  Figure  2. 


Table  1.  Model  Parameters  for  Equation  (4)  by  Fitting  Observed  Potential 
Coefficient  Variation  in  the  Rimge  of  Degree  3  to  36. 


Solution 

A 

B 

RMS 

Residual  x  10^ 

RMS 

%  Residual 

Weighting 

GEM  lOB 

10.6 

2.  17 

0.026 

16.3 

1/obs^ 

GEM  lOB 

13.2 

2.26 

0.014 

16.  8 

equal 

180, 180 

9.41 

2. 12 

0. 037 

16.5 

1/obs^ 

180, 180 

_ 

13.8 

2.30 

0.014 

19.9 

equal 

We  have  also  carried  out  fittings  to  equation  (4)  where  we  have  fixed  B  to 
be  2  or  2.5  exactly  and  solved  for  A  for  the  two  ranges  of  degree  3  to  20  and 
3  to  36.  The  results  arc  given  in  Table  2. 
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Table  2.  Fitting  to  Equation  (4)  Enforcing  Parameter  H 
Using  the  GEM  lOB  Potential  Coefficients 


Range 

A 

B* 

RMS 

Residual  x  10' 

RMS 

%  Residual 

Weighting 

3-20 

7.02 

2.0 

0.086 

22.4 

l/obs' 

3-20 

9.57 

2.0 

0.032 

.50.0 

equal 

3-20 

22.3 

2.5 

0.087 

18.6 

l/obs ‘ 

3-20 

17.8 

2.5 

0.025 

24.0 

equal 

3-36 

6.68 

2.0 

0.070 

19.0 

1 /obs' 

3-36 

9.57 

2.0 

0.024 

53.3 

equal 

3-36 

27.0 

2.5 

0.012 

28.1 

l/obs' 

3-36 

17.8 

2.5 

0.019 

36.6 

equal 

*  Fixed  Value 


If  we  consider  the  degree  range  3  to  20,  the  best  percentage  fit  is  obtained 
when  B  is  2.5.  If  we  consider  the  range  3  to  36,  the  best  ix?recntage  fit  is 
obtained  when  B  is  2.0.  Thus,  the  answer  to  the  question  as  to  how  the  potential 
coefficients  decay  vis  2.0  or  2.5  seems  to  depend  on  the  degree  range  being 
considered.  However,  even  in  the  range  degree  3-20  the  use  of  B  =  2.  5  is  only 
slightly  better  than  the  2.0  value.  To  graphically  display  the  differences,  we 
show  in  Figure  3,  the  case  3  to  36  with  B  =  2  and  B  =2.5  with  weighting  as 
1/obs^.  It  is  clear  that  at  the  higher  degrees,  the  solution  with  B  =  2  gives  :i 
better  fit  to  the  data  than  the  B  =  2. 5  case. 

We  next  consider  the  potential  coefficient  variation  implied  by  the  solution 
to  degree  180.  We  show  in  Figure  4  this  variation  along  with  the  variation  implied 
by  Kaula's  rule  of  thumb  and  the  Tsche ming/Rapp  anomaly  degree  variance  model. 
At  low  degrees,  the  actual  field  has  less  power  than  the  two  models.  The  Tsche r- 
ning/Rapp  model  has  more  power  at  all  degrees,  and  the  Kaula  rule  has  less  power 
after  degree  70.  In  Figure  5,  we  show  again  the  180, 180  variation  along  with  the 
Wagner  and  Colombo  (1978)  model  and  the  five  shell  model  of  Heller  and  Jordan 
(1979).  The  Wagner  and  Colombo  model  consists  (for  the  range  of  considered 
here)  of  five  different  functions  of  the  form  of  equation  (4).  The  Wagner  and 
Colombo  model  fit  our  data  much  bettor  than  the  data  of  Heller  .and  Jordan. 

We  have  fitted  equation  (1)  to  the  data  with  results  shown  in  Table  3  for 
the  case  of  two  parameters  being  found,  and  in  the  case  B  is  fixed  at  2  or  2.  5. 


w  > 


p 


Table  3.  Model  Parameters  for  Equation  (4)  from  Fitting  Observed  Potential 
Coefficient  Variations  in  the  Rjinge  of  Degree  3  to  180 


A 

B 

RMS 

Residual  X 10® 

RMS 

%  Residual 

Weighting 

4.10 

1.80 

0.045 

13.8 

l/obs'’ 

13.7 

2.29 

0.006 

59.  8 

equal 

9.71 

2.0* 

0.011 

26.0 

l/obs‘ 

9.54 

2.0* 

0.011 

26.0 

equal 

81.6 

2.5* 

0.036 

106 

l/obs‘ 

17.7 

2.5* 

0.007 

76.0 

equal 

*  Value  fixed  in  the  adjustment. 


We  show  in  Figure  6  plots  of  the  fits  for  the  first  two  solutions  where  we 
see  that  the  solution  with  the  minimum  percent  residual  fits  best  at  the  higher 
degrees.  In  the  case  of  B  fixed  at  2  there  is  no  substantial  difference  in  the 
results  from  either  weighting  scheme.  When  B  is  Uiken  as  2.5,  the  fits  are 
substantially  poorer  than  when  B  is  2. 


Anomaly  Degree  Variance  Modelling 

Anomaly  degree  variances  can  be  computed  from  potential  coefficients 
using  equation  (3).  (More  rigorous  equations  are  discussed  in  Jekeli  (1978).) 

The  Tscheming/Rapp  model  took  the  form: 

^  ( /  -  2)  ( £  +  B) 

where  B  is  a  positive  integer.  In  addition  to  A  and  B  as  parameters,  an 

additional  term,  s,  was  introduced  such  that  the  sum  of  the  Cjj  times 

would  yield  the  point  anomaly  variance,  Co,  at  the  surface  of  the  earth.  We  have: 

OO 

Co  =  c,  s^^"^  (6) 

In  Tscheming/Rapp,  the  parameters  were  A  =  425.28  mgal^,  B  =  24,  and 
s  =  0.999617.  These  quantities  were  obtained  by  fitting  the  model  (from  degree 
3)  to  observed  anomaly  degree  variances  to  degree  20,  l'’  and  5°  block  anomaly 
variances,  and  a  point  anomaly  variance  of  1795  mgaF'  with  respect  to  an  ellip¬ 
soidal  reference  field. 

For  this  paper,  we  have  extended  this  fitting  to  the  anomaly  degree  variances 
to  degree  180  but  have  excluded  the  block  variances  because  they  have  become 
extraneous  information.  We  have  introduced  information  on  the  horizontal  gradient 
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variance  Gen  for  help  in  strengthening'  the  solution.  The  horizontal  gradient 
variance  is  related  to  the  anomaly  degree  variances  by  (Jekeli,  1978): 


Co.  = 

1  =  2 


Cf  S 


(7) 


(Pertinent  discussions  on  the  interpretation  of  K,  c„,  and  a  may  be  found  in 
Jekeli  (ibid).^  However,  we  have  adopted  for  use  with  equation  (5)  only,  a 
Gqh  =  3500  E  which  is  that  implied  by  the  Tseheming/li'ipp  mt)del.  This  value 
is  too  high  in  terms  of  the  real  world  as  pointed  out  by  Moritz  (1977)  but  as 
shown  by  Jekeli  (1978),  the  Tscheming/llapp  model  is  not  compatible  with  a  low 
horizontal  gradient  variance. 


We  therefore  made  a  parameter  fit  to  the  observed  anomaly  degree  variances 
from  degree  3  to  180,  a  point  anomaly  variance  of  1800 ±25  mgal^  and  a  horizontal 
gradient  variance  of  3500±  50  E^.  The  computer  pmgram  used  was  originally 
written  by  Jekeli  for  his  earlier  studies.  A  number  of  different  runs  wert*  made 
with  different  B  values  and  weighting  schemes.  Tw'o  solutions  of  specific  interest 
are  detailed  in  Table  4.  Besides  giving  Co,  Goh,  Lq  (the  undulation  variance) 
of  the  model,  we  give,  in  the  table,  the  correlation  length  which  is  defined  as 
^  where  C  (O  =  Co/2. 

Table  4.  Parameters,  Residuals,  and  Adjusted  Values  Obtained  in  Kitting 
to  the  Anomaly  Degree  Variance  Model  of  Tscherning/Rapp 


Solution  One 

Solution  Two 

A 

135. 19  mgaP 

429.48  mgar 

B 

5 

24 

s 

0.999780 

0.999613 

RMS  C£  Residual 

1. 35  mgal^ 

3.08  mg!il‘^ 

RMS  %  Residual 

49 

108 

Co* 

838  mgal^ 

1801  mgaf 

Gqh 

3424  e"^ 

3500  E® 

Lo  + 

599  m® 

617  m*" 

C 

65  km 

43  km 

*  Add  7.  6  mgal^  to  refer  to  ellipsoidal  model, 
f  Add  314  m^  to  refer  to  ellipsoidal  model. 


Solution  One  was  made  to  givt'  the  best  overall  fit  to  the  anomaly  degrt'c 
variances.  Solution  Two  was  made  to  fit  the  anomaly  degri'e  variances  as  well  as 
possible  and  to  yield  a  point  anomaly  variance  of  about  1800  mgar.  This  latter 
solution  has  a  B  value  the  same  as  the  original  Tsche ming/Rapp  model,  with  the 
A  and  s  values  being  only  slightly  different.  Because  of  this  agreement,  w'c 
consider  that  the  parameters  used  in  the  earlier  model  are  still  valid  provided 
we  are  willing  to  accept  the  high  gratlient  variance. 
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Those  anoniah  dogn-o  varianoo  models  imply  a  potential  coefficient  variation 
found  by  substituting  (1)  into  (:{).  Plots  of  these  variations  are  shown  in  Figure  7. 
Solution  Two  shows  systematically  more  i>ower  than  the  observed  field  while 
Si)lution  One  shows  less  power  after  about  degree  15. 

In  ortler  to  avoid  the  high  gradient  variance  associated  with  the  model 
given  as  equation  (5)  Moritz  (1977)  suggested  a  two  component  model.  Jekeli 
(197S)  slightly  recast  this  model  into  the  following  form; 

cj  o, - -  a,  +  — — - — — —  oj  8) 

where  Oj ,  CTi,  CT;  ,  A  and  B  are  ciuantities  to  be  determined.  In  practice, 

A  and  B  are  integers  that  are  varied  while  the  other  quantities  are  solved  for  in 
an  adjustment  procedure.  This  procedure  was  implemented  by  Jekeli  (1978)  who 
found  pjirametei-s  based  on  available  data  which  included  the  c.  values  implied 
by  tht'  CKM  9  poU-ntial  cot'fficients,  a  Co.  and  Goh  values.  Tne  c^  values 
given  in  (H)  will  refer  to  a  sphere  whose  radius  is  the  mean  radius  of  the  earth. 
Consequently,  the  jioint  anomaly  variance  on  the  surface  of  this  sphere  can  be 
found  by  substituting  (H)  into  (6)  and  setting  s  to  1. 

The  fitting  to  the  c^  model  given  in  equation  (8)  was  carried  out  using  the 
given  anomaly  degive  variances  to  degree  ISO,  a  Cn  1800  mgals^,  and  a 
C'nn  800  K^.  Initial  fits  with  Gqh  used  by  Jekeli  (1978)  and  Moritz  (1977) 
did  not  yield  go(xl  daUi  fits.  The  value  of  800  was  chosen  considering  the 
observations  described  by  Hein  et.  aU,  (1979).  However,  Goh  still  remains  a 
weak  link  in  our  model  development. 


Approximately  ‘if)  different  runs  were  made  w  ith  different  A  and  B  values 
and  different  wfighting  schemes.  (In  these  fits  the  anomaly  degree  variances 
were  conqxited  from  (J)  w  ithout  any  further  reduction  to  a  mean  sphere.)  We  give 
the  ri'sults  of  these  fits  in  Table  5  for  two  solutions.  Case  One  is  the  solution 
tliat  gives  the  best  overall  fit  to  the  data,  and  Case  IVo  is  the  one  that  best  fits 
the  observed  anomaly  degree  variances.  These  nvo  solutions,  in  terms  of 
potential  coefficient  variations,  are  shown  in  Figure  8.  At  high  degrees,  we  see 
that  Case  One  has  mort'  [lower  than  the  actual  field  or  Case  Two.  This  occurs  in 
this  mixlel  to  satisfy  the  netxl  for  a  Cq  of  about  1800  mgaF. 


FiKurc*  S. 


DEGREE 


Table  5.  Parameters,  Residuals,  and  Adjusted  Values  obtained  in  Fitting 
to  the  Anomaly  Degree  Variance  Model  of  Moritz 


Case  One 

Case  Two 

Oil 

3.4050  mgal^ 

14.966  mgal^ 

Oi  3 

140.03  mgal^ 

999.  25  mgal^ 

.998006 

.987969 

(Ji; 

.914232 

.  850000 

A 

1 

75 

B 

2 

20 

RMS  Cl  Residual 

1.08  mgal'’ 

0.93  mgaP 

RMS  %  Residual 

39.2 

15.  1 

Co* 

1794  mgal"^ 

647  mgaf 

Gqh 

1053  Ef 

15  ir 

Lo  + 

616  m"^ 

641  m^ 

? 

24  km 

110  km 

*  Add  7.6  mgal^  to  refer  to  ellipsoidal  model, 
t  Add  314  to  refer  to  ellipsoidal  model, 
t  f  Constrained  to  this  value. 


Summary  and  Conclusions 

In  this  paper,  we  have  tried  to  improve  our  knowledge  of  potential  coefficient 
and  anomaly  degree  variance  behavior  by  comparing  new  data  to  existing  models 
and  by  obtaining  new  parameters  of  various  models.  The  primary  new  information 
going  into  these  models  is  a  set  of  potential  coefficients  given  to  degree  180  based 
on  satellite  derived  potential  coefficients,  terrestrial  gravity  anomalies,  and 
anomalies  derived  from  Cieos-3  altimeter  data. 

One  of  the  first  steps  was  to  compare  the  potential  coefficient  variations  of  the 
GEM  lOB  solution  to  the  variations  implied  by  the  Kaula  rule  of  thumb  and  the 
Tscheming/Rapp  degree  variance  model.  We  found  that  both  models  have  signifi¬ 
cantly  more  power  than  exists  in  the  real  gravity  field.  Two  parameter  fits  to  the 
data  were  carried  out.  Other  fits  were  carried  out  to  see  if  the  potential  coefficients 
decay  as  l/t®  or  l/f  ■  .  Choosing  tht'  criteria  as  a  best  percentage  fit,  the  best 
fit  is  obtained  with  a  fall  off  of  1//*  for  I  between  3  and  36.  For  I  between  3 
and  20,  a  fall  off  of  \/ is  slightly  better  thjin  l/f‘. 

We  then  made  comparisons  and  fits  to  a  number  of  different  potential  coeffi¬ 
cient  and  anomaly  degree  variance  m«>dels  from  I  equal  3  to  180.  We  summarize 
these  results  in  Table  6  for  a  number  of  common  quantities. 


Table  6.  Summary  of  Model  Comparisons  and  Fits  With  and  To 
the  Potential  Coefficient  Solution  to  Degree  180 


Model 


Number  of 

Pot.  Coeff. 

Variances* 

Correlation 

Parameters 

Residuals 
RMS  X  10®  % 

Co 

(mgal®) 

Gqh 

Lo 

(m^) 

Length 

(km) 

2 

.0126  28 

nil'*’ 

1600+ 

705+ 

3 

.0074  38 

1787 

3500 

612 

42 

10 

.0200  75 

1816 

255 

556 

66 

12 

.0078  14 

1193"^ 

1656+ 

642+ 

2 

.0449  14 

1556^ 

8210+ 

216+ 

2 

.0061  60 

271+ 

28'^ 

594+ 

3 

.0074  39 

1801 

3500 

617 

43 

.0048  25 

838 

3424 

599 

65 

6 

.0072  17 

1794 

1053 

616 

24 

6 

.0640  16 

647 

15 

641 

110 

Kaula 

Tsche  ming/liapp 
Heller/jordan 
Wa  gne  r/C  ol  o  mbo 
A/£®,  case  1  ♦♦ 
A/I®,  Case  2 
One  comp,  f  Case 
C£  model  lease  2^^ 
Two  comp,  r Case  1 
c«  model  lease  2*^ 


+  From  Degree  .8. 

+  ♦  easel;  A  =  4.10x10'®,  B  =  1.80;  Case  2;  A  =  13.  7  x  10~®,  B  =  2. 29. 

ft  Case  1:  A  -  429.48  mgal^,  B  =  24;  Case  2;  A  =  135.19  mgal^  B  =  5. 

H-  See  Table  4. 

f  s  =  0.999017 

Of  the  existing  models,  the  Wagner/Colombo  model  seems  best  although 
the  anomaly  variance  is  less  than  the  1800  mgal^  we  feel  is  reasonable.  Neither 
case  of  the  four  A/ fits  all  quantities  well  although  the  fit  to  potential  coeffi¬ 
cient  variations  is  quite  good  for  case  1.  The  fits  of  the  one  component  degree 
variance  model  do  not  appear  significantly  better  than  the  original  Tscheming/ 
Bapp  model.  The  two  component  model,  case  1,  seems  to  give  the  best  overall 
fit  to  the  daUi  we  have  knowledge  of,  although  the  correlation  length  seems  to 
be  shorter  than  one  usually  expects.  However,  it  agrees  quite  well  with  the 
estimate  of  22  km  obtained  from  the  anomaly  covariance  function  given  in  Table 
19  of  Rapp  (1964). 

We  have  also  computed  a  1°  covariance  function  from  the  two  component 
degree  variance  model  (A  =  1,  B  =  2)  which  is  shown  in  Figure  9.  We  also 
show  the  empirically  determined  l”  covariance  function  from  Tscheming  and 
Rapp  (1974,  Table  15)  with  both  functions  scaled  to  have  the  same  variance. 

The  agreement  is  quite  good,  especially  at  the  location  of  the  zero  crossings. 

We  specifically  conclude  that  for  covariance  studies  not  dealing  with 
accuracy  estimations  in  gradient  computations,  the  parameters  of  the  Tscheming/ 
Rapp  model  are  still  valid.  For  a  better  fit  to  all  expected  values,  the  fit  to  the 
two  component  model  suggested  by  Moritz  (specifically  case  1  in  Table  6)  seems 
to  be  best. 
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Figure  9 


One  Degree  Covariance  Function  From  Empirical  Data  and  as  Derived 
From  the  Two  Component  Anomaly  Degree  Variance  Model  with 

A  =  1,  B  =  2 


SEPflRRTION  (DEGREES) 


30.00 
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